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It is a pleasure for us to dedicate this paper to Prof . Herbert Walther, 
a pioneer in quantum òptics from the very beginning. The investi- 
gation of the amazing properties of single photons both in the mi- 
crowave and the optical domain has always been a central theme in 
his research. We wish him all the best in the years to come! 



I. INTRODUCTION 



Four dècades after the pioneering work on optical coher- 
ence and photon statistics by Glauber (jl, the controlled gen- 
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eration of single photons with well-defined coherence proper- 
ties is now of fundamental interest for many applications in 
quantum information science. First, single photons are an im- 
portant ingredient for quantum cryptography and secure quan- 
tum key distribution |2]. Second, the realization of quantum 
computing with linear òptics (LOQC), which was first pro- 
posed by Knill, Lafiamme and Milburn 1 3], relies on the avail- 
ability of deterministic single-photon sources. And third, var- 
ious schemes have been proposed to entangle and teleport the 
spin of distant àtoms, acting as emitters of single photons, by 
means of correlation measurements performed on the single- 
photon light fields |4j|5l|fi|7j]. Therefore, in recent years, a lot 
of effort has been made to realize single-photon sources. As 
a result, the controlled generation of single photons has been 
demonstrated in various systems, as summarized in a review 
article of Oxborrow and Sinclair |8]. 

Using the process of spontaneous emission from a single 
quantum system is the simplest way to realize a single-photon 
source. In this case, the quantum system is excited by a short 
làser pulse and the subsequent spontaneous decay of the sys- 
tem leads to the emission of only one single photon. This 
has been successfull y de monstrated many times, e.g., using 
sing le molècules Í9Í ll(l \^&, single àtoms fLÏTl . single ions 
11311 . single color cente rs 114 1151(101 or single semiconductor 
quantum dots fl7l fÏ8l fi 9i l20ll - If the quantum system radi- 
ates into a free-space environment, the direction of the emitted 
photon is unknown. This limits the efficiency of the source. 
To overcome this problem, the enhanced spontaneous emis- 
sion into a cavity has been used. The system is coupled to a 
high-finesse cavity and the photon is preferably emitted into 
the cavity mode, which defines the direction of the photons. 
Although a cavity is used, most properties of the photons, like 
the frequency, the duration and the bandwidth, are given by 
the specific quantum system. Only if the generation of single 
photons is driven by an adiabatic passage can these spectral 
parameters be controlled. This technique uses the atom-cavity 
coupling and a làser pulse to perform a vacuum stimulated 
Raman-transition (STIRAP), which leads to the generation of 
one single photon. Up to now, this has been demonstrated 
with single Rubidium àtoms I2jl 12211 . single Caesium àtoms 
1 23] and single Calcium ions 1 24] placed in high-finesse opti- 
cal cavities. 

The characterization of a single-photon source usually 
starts with the investigation of the photons statistics, which 
is done by a correlation measurement using a Hanbury 



Brown and Twiss 1 25 J setup. The observation of antibunch- 
ing shows that the source emits only single photons. How- 
ever, the requirements on a single-photon source for LOQC 
and for the entanglement of two distant àtoms go far beyond 
the simple fact of antibunching. The realization of these pro- 
posals relies on the indistinguishability of the photons, so that 
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even photons from different sources need to be identical with 
respect to their frequency, duration and shape. Therefore it 
is desirable to investigate the spectral and temporal proper- 
ties of single photons emitted from a given source. We em- 
phasize that properties like bandwidth or duration always deal 
with an ensemble of photons and cannot be determined from a 
measurement on just a single photon. Therefore any measure- 
ment of these properties requires a large ensemble of succes- 
sively emitted photons. Several methods have been employed 
to characterize these. 

The first measurement of the duration of single photons has 
been performed by Hong, Ou and Mandel 1 26] . In this ex- 
periment, the fourth-order interference of two photons from 
a parametric down-conversion source was investigated by su- 
perimposing the signal and the idler photon on a 50/50 beam 
splitter. The coincidence rate of photodetections at the two 
output ports of the beam splitter was measured in depen- 
dence of a relative arrival-time delay between the two pho- 
tons. Indistinguishable photons always leave the beam splitter 
together, so that no coincidence counts can be observed. If the 
photons are slightly different, e.g. because they impinge on 
the beam splitter at slightly different times, the coincidence 
rate increases. Therefore, as a function of the photon delay, 
the coincidence rate shows a minimum if the photons impinge 
simultaneously on the beam splitter, and for otherwise identi- 
cal photons the width of this dip is the photon duration. The 
minimum in the coincidence rate goes down to zero if the pho- 
tons are identical. Any difference between the two interfering 
photons reduces the depth of this dip. The first demonstra- 
tion of such a two-photon interference of two independently 
emitted photons from a quantum-dot device has been shown 
by Santoriet al. jz^ . 

In addition to the two-photon coincidence experiments, a 
correlation measurement between the trigger event and the de- 
tection time of the generated photon can be used to determine 
the temporal envelope of the photon ensemble 1 22j |23j, 12411 . 
This latter method is insensitive to the spectral properties of 
the photons. In general, it does not reveal the shape of the 
single-photon wavepackets, unless all photons are identical. 
In case of variations in the photon duration or a jitter in the 
emission time, only the temporal envelope of the photon en- 
semble is observed. No conclusions can be drawn on the en- 
velope of the individual photons. 

The Standard way to determine the coherence time of a 
given light source is the measurement of the second-order in- 
terference using a Mach-Zehnder or Michelson interferome- 
ter. This measurement can also be done with single photons, 
so that each single photon follows both paths of the interfer- 
ometer and interferes with itself. The detection probability 
of the photons at both outputs of the interferometer shows a 
fringe pattern if the length of one arm is varied. The visibility 
of this pattern depends on the length difference of both arms 
and determines the coherence length (or the coherence time) 
of the photons. This method has been used by Santori et al. 
l27ll and Jelezko et al. Í28ll to measure the coherence time of 
their single-photon sources. However, this method is hardly 
feasible for photons of long duration, because the length of 
one arm of the interferometer must be varied over large dis- 



tances. Furthermore, the measurement depends on the me- 
chanical stability of the whole setup, i.e. the interferometer 
must be stable within a few per cent of the wavelength of the 
photons, which might not be easy. 

Only recently, adiabatic passage techniques have allowed 
the generation of photons which are very long compared to 
the detector time resolution. Therefore the detection time of 
a photon can be measured within the duration of the single- 
photon wavepacket. As a consequence, the two-photon inter- 
ference can be investigated in a time-resolved manner, i.e., 
the coincidence rate can be measured as a function of the 
time between photodetections 1 29, 30]. The theoretical anal- 
ysis shows that this method not only gives information about 
the duration of single photons, but also about their coherence 
time. Here we discuss how to use this method for a spectral 
or temporal characterization of a single-photon source. 

The article is organized as follows: After a brief summary 
of the nature of single-photon light fields (Chapter|HJi, we dis- 
cuss the interference of two photons on a beam splitter and 
introduce the time-resolved two-photon interference (Chapter 
llllt . Thereafter, we show how a frequency and an emission- 
time jitter affects the results of a time-resolved two-photon 
interference experiment (Chapter llV> . On this basis, the ex- 
perimental characterization of a single-photon source, based 
on an adiabatic passage technique, is discussed (ChapterlvV 

II. SINGLE-PHOTON LIGHT FIELDS 

The quantum theoretical description of light within an op- 
tical cavity is well understood 1 3 1 ] . The electromagnètic field 
between the two mirrors is subject to boundary conditions 
which lead to a discrete mode structure of the field. Each 
mode can be labeled by a number l and is characterized by 
its individual frequency, lo. These eigenfrequencies are sepa- 
rated by Aw = 2ttc/L, where L is the round-trip length of the 
cavity. The quantization results in a discrete set of energies, 
E n = huj(n+l/2) and the appropriate eigenstates are defined 
by means of creation, al, and annihilation, à t , operators. The 
energy eigenstates 



are states with a fixed photon number, n. In this context, pho- 
tons are the quanta of energy in the modes of the cavity. 

In the limit of L — > 00 and Aw — > 0, the mode spectrum 
becomes continuous. In this case it is convenient to introduce 
continuous-mode operators à^(u) and à(u>) according to 

a\ -»• (AwJ^óHu) (2) 

These operators create and annihilate photons as quanta of 
monochromatic waves in free space. These waves of infinite 
spatial extension do not have any beginning or any end. How- 
ever, photons generated in a laboratory are characterized by 
a certain frequency bandwidth or a finite spatial extension. 
Therefore it is desirable to define operators which create or 
annihilate photons in modes of a given bandwidth, or, in other 
words, of a well defined spatiotemporal spread. 
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A. Frequency modes 

In contrast to modes describing monochromatic waves, it is 
possible to define field modes of a given frequency distribu- 
tion. These modes represent wavepackets travelling with the 
speed of light c through the vacuum, and the bandwidth k of 
such a mode determines the duration St of the wavepacket. A 
frequency distribution is described by a normalized complex 
function which is called the mode function of the field. 
The operators à'(u>) and à(u>) can be used to define a new 
set of operators for the creation and annihilation of photons in 
these new modes l32ll . The creation operator, e.g., is given by 



duj x(w) ér (w). 



(3) 



Note that the mode function x( w ) can be written as the 
product of a real amplitude, e{ui), and a complex phase, 
exp (— The phase term includes the emission time 
To and the propagation of the wavepacket. In the following, 
we restrict our discussion to Gaussian wavepackets centered 
at the frequency Uq. Their mode functions read 



(üJ - UJ f 



exp exp(-icj(r + z/c)). 

\ K J 

(4) 

For an ideal single-photon source which always produces 
identical photons, the light field is always described by the 
same quantum mechanical state vector. In other words, the 
state vector is given by the creation operator w (x) acting on 
the vacuum state |0) for every single photon: 



1, 



K io>. 



(5) 



We emphasize that such an ideal source is hardly feasible. 
Usually the generation process cannot be controlled perfectly 
and therefore the mode function is subject to small variations. 
To take this into account the light field must be described by a 
quantum mechanical density operator 



g= I \l xW )(l xW \ 



(6) 



Here we assume that the source produces single photons with 
a Gaussian frequency distribution and the parameters of this 
distribution are subject to small variations, according to a dis- 
tribution function /(#). The parameter i? stands for the center 
frequency, ujq, the bandwidth, k, or the emission time, tq, of 
the photon, or a combination of these. 



B. Spatiotemporal modes 

Due to the Fourier theorem, each mode with a certain 
frequency distribution x( w ) can t> e assigned to a temporal 
wavepacket which is travelling through space. A mode with 
the Gaussian frequency distribution given by Eq. @ there- 
fore belongs to a spatiotemporal mode £(f — z/c) of Gaussian 



shape. With the substitution q :- 
by the function 



t — z/c this mode is given 



= \fïà? CX P (~^) exp ( iuJ o{ra - q)) 



e(q) exp(iw (r - q)). 



(7) 



The duration St of this Gaussian wavepacket is given by the 
reciprocal bandwidth of the frequency distribution, St = 2/n. 
Blow et al. 1 32] have shown that creation and annihilation op- 
erators can also be assigned to spatiotemporal modes. In order 
to do that, one has to define the Fourier-transformed operators 



à\q) = (2tt) 



-1/2 



à(q) = (2tt)-^ 2 / dwà(u)e iwq 



(8) 



(9) 



By means of these operators we define a flux operator 
à'(q)à(q). Its expectation value has the unit of photons per 
unit time. We need this operator in the next section to describe 
the detection of single photons. 

Eqs. l|8} and l|9} are only vàlid if the bandwidth of the modes 
is much smaller than the frequency of the light, k << ujq. 
This also limits the localization of a single photon in such 
a spatiotemporal mode. In case of optical freqüències, this 
condition is usually fulfilled. 

In analogy to Eq. the Fourier-transformed operators can 
be used to define creation and annihilation operators for pho- 
tons of spatiotemporal modes £(<?): 



\= / dq£{q)àHq) 



(10) 



To take fluctuations into account, one can again write the den- 
sity operator of the light field as in Eq. 0, but using spa- 
tiotemporal modes. In this case, $ stands for any combination 
of uj , St, and t . 



C. Single-photon detection 

Choosing spatiotemporal modes for describing the state of 
a single-photon light field simplifies the formal description of 
the detection of a photon. We assume a detector with quantum 
efficiency 77 placed at the position z = 0. The response of the 
detector within a time interval [to — dto /2, íq + dío /2] is given 
by the expectation value of the flux operator: 

/■to -dt /2 

P {1) {t )=r] / dttr[Qà\t)à[t)] (11) 

Jto+dto/2 

In case of single-photon wavepackets, the function P^(to) 
gives the probability to detect this photon within the consid- 
ered time interval. In practice, the lower limit of the duration 
dío is given by the detector time resolution T, i.e., dto > T. 

If the photon duration is much longer than the detector time 
resolution, St » T and dío = T, Eq. dl 1> can be simplified 
to 



P (1) (í ) = vTtr[g^(to)à(t ) 



(12) 
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The measurement of the detection probability requires a large 
ensemble of single photons. In the following, we therefore 
assume a periòdic stream of single photons emitted one-after- 
the-other, so that the photons always hit the detector one by 
one. If all photons of this stream are identical, the light field 
can simply be described by a state vector 1 1^) and the density 
operat or is given by g = |lç)(lç|, with = c||0). In 
this case, the average detection probability of the ensemble 
of photons is given by the square of the absolute value of the 
mode function, £(q), and is therefore identical to the shape of 
each individual photonic wavepacket, 

P (1 \t Q )=vT\t;(to)\ 2 =vTe 2 (t Q ). (13) 

As already discussed, the photons may differ from one an- 
other, and the density operator is given according to Eq. 0. 
The average detection probability is then given by 

pW(to)=r)T J dï?/(ï?)e 2 (í ,i?). (14) 

To obtain this equation, we assume that trace and integration 
can be exchanged. Obviously the average detection probabil- 
ity for the photon ensemble differs from that for individual 
photons. The average detection probability is, in general, af- 
fected by the variation, /($), of the parameters of the mode 
function, Therefore it shows only a temporal envelope 
of the photon ensemble. However, the effect of each param- 
eter onto P^(to) can be very different. A variation of the 
frequency, e.g., does not affect the real amplitude of the mode 
function, so that the average detection probability, P^'fto), 
is simply given by Eq. (1131 . This is not the case for variations 
of the other parameters, as will be shown in Chanter lIVl 

III. TWO-PHOTON INTERFERENCE 

We now consider two independent streams of Gaussian- 
shaped single photons that impinge on a 50/50 beam split- 
ter such that always two photons are superimposed. As we 
show in Chapter|V] these two streams can originate from one 
single-photon source by directing each photon randomly into 
two different paths of suitable length, so that successively gen- 
erated photons hit the beam splitter at the same time. Here we 
ask for the probability to detect the photons of each pair in 
different output ports of the beam splitter. In case of iden- 
tical photons, the joint detection probability is zero. With 
polarization-entangled photon pairs emitted from a down- 
conversion source, this effect has first been used by Alley and 
Shih 13311 to test the violation of Bell's inequality by joint pho- 
todetections, and one year later, Hong, Ou and Mandel lEóTl 
have used it to measure the delay between two photons with 
sub-picosecond precision. Recently, two-photon interference 
phenomenona have successfully been employed to test the in- 
distinguishability of independently generated single photons 
ll27Íl . To illustrate this interference effect, we first assume that 
each photon of a given stream can be described by the same 
quantum mechanical state vector, 1 1^), but allow the state vec- 
tors of the two considered streams to differ from one another. 



In Chapter lTVl we generalize this discussion to streams of pho- 
tons which show a variation in the parameters of the mode 
functions, e.g. a variation in the photon frequency. Finally we 
show that the interference of photon pairs reveals information 
about the variations of the mode functions. 

In section llll Al and lTll Bl we start with a brief discussion of 
the beam splitter and the principle of the two-photon interfer- 
ence. Afterwards we analyze the joint detection probability 
for photons in the limits of a photon that is either very short 
or very long compared to the detector time resolution. 



A. Quantum description of the beam splitter 

The beam splitter is an optical four-port device with two 
inputs and two outputs. The principle of the beam splitter 
is shown in Figure[2 As discussed by Leonhardt [34], each 
port has its own creation and annihilation operators, and the 
output operators can be expressed by the input operators using 
a unitary transformation matrix B. This relation is vàlid for 
creation and annihilation operators à(u) of monochromatic 
waves as well as for operators of spatiotemporal modes, b x or 
cj. Itreads: 

ft)= B ft) ^ («UH«U)B* (15) 

In the following discussion, we assume an ideal lossless and 
polarization independent beam splitter with transmission co- 
efficient ^fa. The matrix of this beam splitter is given by 

The opposite signs of the off-diagonal terms reflect the phase 
jump of 7r for the reflection at one side of the beam splitter. 

The transmission of photons from the input side to the out- 
put side of the beam splitter can be understood as a quan- 
tum mechanical evolution of the system. This evolution can 
be described in two equivalent pictures corresponding to the 
Heisenberg and the Schròdinger picture in quantum mechan- 
ics Í35l 13611 . In the Heisenberg picture, the evolution is de- 
scribed by the creation and annihilation operators. The out- 
put operators are considered as the evolved input operators 
whereas the state vector of the field remains unchanged (see 
Figure Q ( a ))- Using the unitary Operators B and B\ this 
evolution can also be expressed by 




Alternatively, in the Schròdinger picture, the evolution can be 
calculated using the state vector of the light field. In this case, 
the state vector of the input side | evolves to a state vector 
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(ü) Heisenberg 
1 




(b) Schródinger 
1 



fit 



FIG. 1: The ideal lossless beam splitter is fully characterized by its 
transmission coefficient a. The reflectíem coefficient is then given 
by \/l — a. Light can enter the beam splitter through two different 
input ports 1 and 2. According to the transmission and the reflection 
coefficient, it is divided into the output ports 3 and 4. For one of the 
reflections, the light is subject to a phase jump of n which is indicated 
by the minus sign. In the Heisenberg picture (a), one accounts for 
this process by transforming the creation and annihilation operators 
of the two input modes (1 and 2) into suitable operators of the output 
modes (3 and 4), whereas in the Schródinger picture (b), the action 
of the beamsplitter is expressed by the unitary operator B\ which 
acts on the wavevector and couples the two through-going modes ( 1 
and 2). 



(a) ! 

2 



+ 



(b) , 

2 



(c) , 

2 



(d) ! 

2 



+ 



FIG. 2: Two impinging photons lead to four possible photon distri- 
butions at the beam-splitter output. In the first two cases (a) and (b) 
the photons would be found together. In the remaining two cases (c) 
and (d) the photons would leave the beam splitter through different 
ports. Since the quantum states of the cases (c) and (d) show different 
signs, they interfere destructively. 



at the output side, |^ uí) = B^\^i n ), while the modes them- 
selves do not change, that is modes 1 and 2 are defined as the 
through-going modes (see Figure[2(b)). In the next section 
the Schródinger picture is used to illustrate the principle of 
the two-photon interference. 



B. Principle of the two-photon interference 

We consider two identical photons that impinge on a 50/50 
beam splitter. The input state of the light field is given by 
l^m) = |l)2ír, where the indices label the two input 

ports and the polarization of the photons. Here, we assume 
two photons of horizontal polarization. In the Schródinger 
picture we describe the evolution of the state using the unitary 
operator B' as follows: 



B^\í) 1H \l) 2H 



D U 1H U 2H 



Since BB^ is equal to the identity operator 1 and B^ |0) = |0) 
we can write 



B ] à\ H à\ H \Q) = &à\ H B &à\ H B\Q), 
which according to Eq. (II Vi gives 



B^à\ H B B^à\ H B\Q) 

= VV2 (à\ H - à\ H ) 1/V2 (à\ H + à\ H ) 



1/2 (à\ 2 H - àH + 



|0) 
|0>. 



(19) 



(20) 



Each term in this sum of creation operators corresponds to 
one of four possible photon distributions in the beam splitter 
output ports, shown in Figure[2] In the first two cases, (a) and 
(b), both photons are found in either one or the other output, 
whereas in the cases (c) and (d), the photons go to different 
ports. The last two cases are indistinguishable, but the two 
expressions leading to cases (c) and (d) have opposite sign. 
Therefore the two possibilities interfere destructively. As a 
consequence, the two photons always leave the beam splitter 
as a pair and the output state is given by the superposition 



St|l) lff |l) 2H = -L(|2) m |0) 2ff - \0) 1H \2) 2H ). 



(21) 



This quantum interference oceurs only if the photons are iden- 
tical. If the photons were distinguishable, no interference 
takes place. For example, two photons of orthogonal polar- 
ization, l^m) = |1)iít|1)2V, give rise to four different output 
states which are distinguishable by the photon polarization. 
In this case the overall output state can be written as a product 
state, e.g., 



fí t |l)m|l)2v = 
^(|l)m|0) 2 -|0) 1 |l) 2íí )® 



(|l)iv|0) 2 + |0)i|l) 2 y), 



which describes the state of two independently distributed 
photons. 

Note that all temporal aspeets of the light field are neglected 
in the above discussion. In the next chapters, the two-photon 
interference is discussed under consideration of the photon 
duration and the time resolution of the detectors. 



C. Temporal aspeets of the two-photon interference 

We now take into account that the photodetections in the 
output ports of the beam splitter might oceur at different times, 
t\ and t 2 . We use the Heisenberg picture to calculate the prob- 
ability of a joint photodetection from the second-order corre- 
lation function, 



G (2) (íi,í 2 ) =^]tr[à,2 A 3sAs ,(h,t 2 ) }, 



(22) 



where gi i2 describes the two-photon input state and the opev- 
ator A 3sAs ,(t 1 ,t 2 ) is given by 



|0) 



(18) 



À 3sAs ,(t 1 ,t 2 ) := 4 s ( í i)«L'( í 2)a 4s '( í 2)ò 3;5 ( í i) 
and s,s'e{H,V}. 



(23) 
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The probability for a photodetection at the first detector within 
the time interval [í — dto/2, t + dío/2] and at the second 
detector within a time interval shifted by r, [ío+r— dr/2, ío+ 
t + dr/2], is then given in analogy to Eq. dl Q , 



to+dto/2 ío+r+dr/2 



P (2) (W 



J dh J dí 2 G^Ht^h). (24) 



ío-dto/2 ío+r-dr/2 



Here we assume that the detectors have different efficiencies, 
?73 and 774. In analogy to Section lïTcl the smallest duration of 
the detection intervals is given by the detector time resolution, 
T, so that dío > T and dr > T. In the following, we calculate 
the joint detection probability in the limit of very short and 
very long photons. 

If the photons are very short compared to the time resolu- 
tion of the detectors, St «T, the limits of the integration in 
Eq. i24\ can be extended to infinity, so that 

P (2) =M4 JJdhdhG^ituh) (25) 

gives the probability of a coincidence of photodetections. 

For very long photons with St >> T, the integration in 
Eq. Í241 leads to 



P (2 Ht Q ,T) = 773774 G (2) (*o.*o +t) dí dr 



(26) 



Therefore the probability of a joint photodetection can be 
studied as a function of the two detection times, íq an d ío + t. 
In practice, only the time difference, r, between two photode- 
tections is relevant. Therefore we integrate P( 2 '(ío,r) over 
the time ío of the first photodetection. This gives 



P {2) (r)=r )3 r H T / dí G^(í , h + r), 



(27) 



where dr is substituted by the detector time resolution T. The 
second-order correlation function, G^ 2 \ plays a central role in 
the calculation of the joint detection probability. It is now an- 
alyzed taking the polarization and the spatiotemporal modes 
of the photons into account. 



cos<^ |lç 2 )2Jï|0)2y + sinip |0) 2 íï|1ç 2 
density operator, g 12 = j*i n )(*w»l' 
state 



2V, respectively. The 
is given by the input 



cosíp m i ) 1H \í (2 ) 2 H + simp \1£ 1 )ih\M 2 )2V, (28) 



which is a superposition of the cases in which the impinging 
photons are parallel and perpendicular polarized to each other. 
The correlation function can then be written as a sum of two 

^,(2) , „(2) 

expressions G HH and G HV , 



G {2) = cos 2 ip G 



(2) 

HH 



sin 2 (p G (2) 



(29) 



(2) 

where the first function G HH accounts for the input state in 
which both photons have the same polarization and the second 
function G HV accounts for the perpendicular polarized case. 
Taking the mode functions into account, these two expressions 
read 



G% ) H {t 1 ,t 2 ) 



I6(íi)£2(í2)-6(íi)£i(t 2 )| 2 



(30) 



r (2) ,. . v ia(^i)6ft2)i 2 + iafe)6fti)i 2 nn 



We emphasize that the correlation function for parallel polar- 
ized photons is always zero for t\ = t 2 , even if the mode 
functions £1 (í) and £2^) are not identical. As a consequence, 
the probability of a joint photodetection, Eq. J27L is always 
zero for r = Í2 — íi = , i.e. no simultaneous photodetec- 
tions are expected even if the photons are distinguishable with 
respect to their mode functions. 

As already mentioned in Cbapter flI Al the mode function 
can be written as the product of a real amplitude and a com- 
plex phase, £j(í) = £j(t) exp (— i $j(í)) with j 6 1,2. 

(2) 

Since the correlation function G HV for perpendicular polar- 
ized photons is independent of the phase, it can be written as 



D. Correlation function 

We calculate the correlation function G^ for two photons 
characterized by two mode functions, £1 and £2- Without loss 
of generality, we assume that both photons are linearly po- 
larized with an angle <p between the two polarization direc- 
tions. The state of the photons is then given by \1^)ih and 



G^y (íl, Í2) 



(ei(íi)e 2 (Í2)) 2 + (ei(Í2)e 2 (íi)) 2 



(32) 



This is not the case for the correlation function of parallel po- 
larized photons which cames a phase-dependent interference 



term. It can be expressed as gSjj(íi,Í2) = G K ^ ! v {t\,t2) 
F(íi,í 2 ), with 



.(2) 



F[tiM) := e^mt^mt,) cos ($i(íi) _ $i(ía) + $2(í2) _ $2(íi)) _ 

1 



However, this phase-dependency is only relevant, if $1 (í) and $2 (t) display a different time evolution. Otherwise the sum 
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over the phases is always zero. Such a difference in the time 
evolution is given if, e.g., the freqüències of the photons are 
different. In that case, the interference term oscillates with 
the frequency difference, which gives rise to an oscillation in 
the joint photodetection probability, P^ (r). This will be dis- 
cussed further in Section lTlI FI 

Taking the interference term into account, the overall cor- 
relation function, Eq. \29\ . can be summarized to 

G( 2 )(Í!,Í 2 ) = G { g> v {t l ,t 2 ) - cos 2 p t 2 ), (33) 

where the effect of the interference term depends on the angle 
ip between the two photon polarizations. 

In the next two sections, the joint detection probability, 
Eq. ( 1241 . is analyzed for very long and very short photons. 



E. Two-photon interference without time resolution 

First we assume Gaussian-shaped photons which are very 
short compared to the time resolution of the photodetectors, 
St << T. In this case one can only decide whether there is 
a coincidence of detections within the time interval T or not, 
and the coincidence probability is given by Eq. (1251 . With a 
possible frequency difference A := QJ02 — <^oi and an arrival- 
time delay St := tq2 — Tqi of the photons, the coincidence 
probability is given by 



P 



(2) 



1 — cos <p exp 



,)l 2 



4/A 2 



exp 



' St 2 



(34) 

where we assume that the photons hit perfect photodetectors 
with 773 = 774 = 1. We analyze the coincidence probability as 
a function of the photon delay St for different photon polar- 
izations and frequency differences. This is shown in Figure|3] 

As already discussed in Section Ull Bl perpendicular po- 
larized photon pairs, ip = n/2, show no interference at all. 
Therefore, the probability for detecting photons at different 
output ports of the beam splitter is always 1/2, independent 
of the photon delay, St. 

If the photons have identical polarizations, ip = 0, and iden- 
tical freqüències, A = 0, the coincidence probability shows 
a Gaussian-shaped dip centered at St = 0. The minimum 
of P( 2 ) is zero, indicating that the photons never leave the 
beam splitter through different output ports. If the photon 
pairs show any difference in their polarization, shape or fre- 
quency, there is no perfect interference and the minimum of 
the dip is no longer zero. Therefore the two-photon interfer- 
ence can be used to test the indistinguishability of photons. 

The first measurement of the coincidence rate as a function 
of the relative photon delay was performed by Hong, Ou and 
Mandel 1 26] using photon pairs from a parametric downcon- 
version source. They controlled the relative delay of the pho- 
tons by shifting the position of the 50/50 beam splitter. The 
freqüències and bandwidths of the photons were adjusted by 



using two identical optical filters, so that the coincidence rate 
dropped nearly to zero for simultaneously impinging photons. 
As one can see from Eq. j34l >. the width of the dip is identi- 
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FIG. 3: Coincidence probability as a function of the arrival time de- 
lay, St, of two linear polarized photons. In case of perpendicular 
polarized photons (tp = n/2) there is no interference at all and the 
coincidence probability shows the constant value 1 /2. If the photons 
are parallel polarized (tp = 0) and have identical frequency (A = 0), 
there is a Gaussian-shaped dip which goes down to zero for simulta- 
neous impinging photons, St — 0. Any difference in polarization or 
frequency leads to a reduced depth of this dip. 



cal to the photon duration, St. Therefore this experiment was 
used to measure the duration and bandwidth of the photons. 

So far, most two-photon interference experiments were per- 
formed with very short photons. Therefore the joint detection 
probability was only considered as a function of the photon 
delay, St. However, if the photon duration is much larger 
than the detector time resolution, the time r between the pho- 
todetections in the two output ports can be measured and the 
joint detection probability can additionally be analyzed in de- 
pendence of the detection-time difference. 



F. Time-resolved two-photon interference 

We now assume, that the photon duration is much larger 
than the detection time resolution, St >> T. Since the time, 
t, between the photodetections can be measured within the 
photon duration, the joint detection probability can be ana- 
lyzed as a function of this detection-time difference. Using 
Eq. 1271 and assuming Gaussian-shaped photons of identical 
duration, St, the joint detection probability is given by 
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(a) no interference: 

perpendicular polarization 



(b) interference: 

parallel polarization 



(c) interference of photons 
with different freqüències 

<0> 




FIG. 4: Joint detection probability, P' 2 ^ , as a function of the relative delay between the photons, St, and the time difference between photode- 
tections, t, for perpendicular polarized photons (a) and parallel polarized photons (b). In Figure (c) the parallel polarized photon pairs have a 
frequency difference A, which leads to an oscillation in the joint detection probability. All times and freqüències are normalized by the photon 
duration, St. 



P&(t,St) = 



T 



/ttíí 



1 — cos 2 (p cos (A r) 



sinh 



r 



st 2 



■ exp 



St 2 



St 2 



(35) 



In Figure@]the joint detection probability is shown as a func- 
tion of the photon arrival-time delay, St, and the detection 
time difference, r. 

The sign of r indicates which detector clicks first. Simi- 
larly, the sign of St determines which photon arrives first at 
the beam splitter. Note that the joint detection probability can 
only be different from zero if |t| sa |<5r|. This leads to the 
cross-like structure in Figure|4](a-c). Since the photons only 
interfere if the relative delay is smaller than the photon dura- 
tion, we focus our attention to the center of Figure@](a-c). 

Again, we start our analysis with perpendicular polarized 
photon pairs. Obviously, no interference takes place, and 
as one can see in Figure|4](a), even simultaneously imping- 
ing photons (with St = 0) can be detected in different out- 
put ports of the beam splitter. The joint detection probabil- 
ity shows therefore a Gaussian-shaped peak. According to 
Eq. i35\ . the width of this peak is identical to the photon du- 
ration, St. Since the photons are distinguishable by their po- 
larization, an additional frequency difference, A, does not af- 
fect this result. Assuming photon pairs with identical mode 
functions, the joint detection probability of perpendicular po- 
larized photons can be used to determine the photon duration. 

Figure|4](b) shows the joint detection probability for par- 
allel polarized photons of identical frequency, A = 0. For 
simultaneously impinging photons the joint detection proba- 
bility is always zero, which indicates that the photons coalesce 
and leave the beam splitter always together. 



If the parallel polarized photons show a frequency differ- 
ence, the joint detection probability oscillates as a function of 
the detection time difference, r. This is shown in Figure|4](c). 
As one can see from Eq. i35\ . the frequency difference, A, 
determines the periodicity of this oscillation. We emphasize 
that the oscillation always leads to a minimum at r = 0, inde- 
pendent of A, so that even photons of different freqüències 
are never detected simultaneously in different output ports. 
Furthermore, the joint detection probability at the màxima is 
always larger than the joint detection probability for perpen- 
dicular polarized photons. 



Without time resolution, the detection-time difference 
cannot be measured and the joint detection probability, 
p( 2 ' (t, St), has to be integrated over r. This links the results 
of a time-resolved two-photon interference to the discussion 
of Section IÏII El In case of perpendicular polarized photons, 
the r-integrated function P^ 1 (St) shows the constant value 
1/2. If the photons are identical, the integration leads to a 
Gaussian-shaped dip, which was already discussed in Section 
IIII El However, the oscillation in the joint detection probabil- 
ity for photon pairs with a frequency difference is no longer 
visible. The integration leads, in accordance to Eq. ( I34l i. only 
to a reduced depth of the dip in P^ (St). 
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IV. JITTER 



A. Frequency jitter 



Up to now, we assumed that all photons of a given stream 
can be described by the same state vector However, this 
requires a perfect single-photon source, which is able to gen- 
erate a stream of photons without any variation in the param- 
eters of the Gaussian mode functions. Here, we consider a 
more realistic scenario, in which a stream of single photons 
shows a jitter in the parameters, #. The quantum mechanical 
state of the photons is then given by the density operator of 
Eq.©. ' 

Such a jitter in the mode functions has important conse- 
quences on the results of measurements which can be per- 
formed on the single-photon stream. On the one hand, as al- 
ready discussed in Section lïïCl it affects the average detection 
probability of the photons in a way that its measurement does 
in general not reveal information about the duration or shape 
of each single photon. On the other hand, variations in the 
mode functions of photon pairs have an infiuence on the joint 
detection probability in two-photon interference experiments. 
This is discussed in some detail in the following two sections. 

To analyze the effect of jitters on the two-photon interfer- 
ence, we consider two streams of Gaussian-shaped photons 
with a variation in the parameters of their mode functions. In 
analogy to Eq. (|6j the density operator for photon pairs im- 
pinging on the beam splitter is given by 



Ql,2 



J f d*! M 2 hi^) / 2 (tf 2 ) |l ei )|l í2 )(l ei |<l e2 |, 



(36) 



so that, using Eq. d22l . the correlation function reads 

G^(t ,t +T) = 

JJd^Mi /i(0i)/ 2 (0 a ) tr(è(£i,£ 2 ) À(t ,t +r)). 

(37) 

Here the expression éíCiíò) substitutes jl^) |lç 2 )(l{ 1 1 (lf 2 1. 
In Eq. (I36> we assumed that all photons have identical polar- 
ization and that they are completely independent from each 
other. Therefore the density operator has only diagonal ele- 
ments. The parameters of the mode functions £i and £ 2 of the 
two streams are summarized by í?i and $2, respectively. In 
general, all parameters of the mode functions could be subject 
to a variation, and all the variations could in principle depend 
on each other. However, in the following two sections, we fo- 
cus our attention only on two examples of jitters and analyze 
the detection probability of photons, P' 1 ', for a single photon 
stream as well as the joint detection probability, P^ 2 \ of two 
streams superimposed on a beam splitter. 

First, in Section IIV Al we consider streams of photons 
which are characterized by a variation of the center frequency, 
ojQj, so that each photon pair exhibits a variation of the fre- 
quency difference, A = 0^02 — All remaining parameters 
of the mode functions, e.g. the duration of the photons, are as- 
sumed to be identical. 



In Section lrV Bl we consider photons, 
which show only a variation in their emission time, so that 
photon pairs are characterized by a variation in their arrival- 
time delay Sr = T02 — tqi . 



We start our discussion of a frequency jitter by analyzing 
its effect on the average detection probability, P^(to), for 
a perfect photodetector with the detection efficiency 7/ = 1 . 
If the frequency variation in the stream of single photons is 
described by a normalized distribution function, f(uj), the av- 
erage detection probability is given, according to Eq. ( H4> . by 
the integral 



P (1) (*o)=T J dw/H |í(í ,^)| 2 



(38) 



Since only the phase of the Gaussian mode functions depends 
on the frequency, the absolute value, |£(ío, i^)| 2 = e 2 (ío)> is 
independent of lo. Thus, the average detection probability is 
not affected by any frequency jitter and is entirely determined 
by the spatiotemporal mode function of each single photon. 

However, a frequency jitter affects the joint detection prob- 
ability of photon pairs superimposed on a beam splitter. To 
illustrate this, we assume two independent streams of pho- 
tons, each fiuctuating around a common center frequency ujq 
according to a normalized Gaussian frequency distributions, 
/i(woi) and f 2(^02)- Hence, the frequency difference of the 
photon pairs, A = UJ02 — ujqi, shows also a normalized Gaus- 
sian variation, 



/(A) 



1 



/7F Suj 



exp(— A 2 /Slü' 2 



(39) 



with width Suj depending on the widths of the frequency dis- 
tributions of both streams, Sui = yj Sui^ + Slüq 2 - The density 
operator of the photon pairs can then be expressed in terms of 
the distribution function of the frequency difference, 



éi, 2 = í dA/(A) g(a,6) 



(40) 



As the operations of trace and integration are exchangeable, 
the correlation function, according to Eq. Í311 . can be written 



G^(t , T ) = / dA/(A) tr(è(6,6)i(ío,ío + r)). (41) 



For photons which are very long compared to the detec- 
tor time-resolution, the joint detection probability is given 
by Eq. i27\ . In case of simultaneously impinging photons, 
St = 0, this leads to 



P( 2 )(r) 



2Vtt<5í 



1 — cos 2 ip exp 



(42) 



For photons of parallel polarization, the result is shown in Fig- 
ure[5] In the limit of Suj — > 00 the joint detection probability 
shows a Gaussian-shaped peak of width T\ = St, which is 
the photon duration. As one can see from Eq. Í421 . the joint 
detection probability is always zero for r = as long as the 
width of the frequency distribution, Suj, is finite. In fact, as 
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FIG. 5: Joint detection probability as a function of the detection- 
time difference, r, for simultaneously impinging photons, St — 0, 
of identical polarization. The photons are subject to a frequency jitter 
of width 8lü. 



one can deduce from Eq. (142b . this leads to a dip in the joint 
detection probability around r = that is Ti = 2/5u> wide. 
Note that represents a coherence time which must not be 
mixed up with the duration of each single photon, St. 

To determine the amount of a frequency jitter from a time- 
resolved two-photon interference experiment, one has to per- 
forin two measurements. First, the joint detection probabil- 
ity of perpendicular polarized photons, ip = n/2, reveals the 
photon-duration, St. Afterwards, the joint detection probabil- 
ity of parallel polarized photons, ip = 0, is used to measure T2 
and derive the width of the frequency jitter, Seu. This is shown 
in detail in ChapterlVl 

If the photons are very short compared to the detector time- 
resolution, the coincidence probability must be calculated us- 
ing Eq. fl25t . The coincidence probability is then a function of 
the relative photon delay, St, and is given by 



<o\ , - •> 1 / 2 cos 2 tp 

P {2) (St) = - 1 - = exp 



TTo • (43) 



St 2 
' St 2 



In analogy to Sec. IIH El a frequency jitter, Su>, now leads to a 
decreased depth of the Gaussian-shaped dip, while the width 
of this dip is not affected and always identical to the photon 
duration. 

In principle, it is possible to derive the frequency jitter 
also from a two-photon interference experiment without time- 
resolution, but there are some major disadvantages. First, the 
depth of the dip depends not only on a frequency jitter, but 
also on the mode matching of the transversal modes of both 
beams. A nonperfect mode matching leads to a factor com- 
parable to cos 2 ip in Eq. ( I43l l. Therefore, in contrast to the 
time-resolved measurement, one cannot distinguish between 



a nonperfect mode matching or a frequency jitter. Second, in 
case of two independent streams of photons from two differ- 
ent single-photon sources, it is impossible to decide whether a 
constant frequency difference or a frequency jitter is the rea- 
son for a decreased dip depth. And third, if the frequency jitter 
is large, the depth of the dip is very small, whereas in a time- 
resolved measurement, the dip-depth remains unchanged. As 
it is much more reliable to determine a small width rather than 
a small depth, the time-resolved method is much more power- 
ful. 



B. Emission-time jitter 

Now we assume a stream of single photons which shows 
a jitter in the emission time of each photon. This variation 
of the emission time is assumed to be given by a normalized 
Gaussian distribution function, /(to). The average detection 
probability of the photons for an ideal photodetector with r\ — 
1 is again given by Eq. H41 . 



PW(Í ) = T |dr/(r) 



(44) 



This is a convolution of the detection probability, |£(ío)| 2 > of 
each single photon and the emission-time distribution, J(tq), 
of the photon stream. Therefore the average detection proba- 
bility is always broader than the detection probability which 
would arise solely from the duration of each single photon. 
This shows that a variation in the parameters of the spatiotem- 
poral mode functions can alter the detection probability of the 
photons. Therefore, in general, the average detection proba- 
bility is not identical to the detection probability of individual 
photons. 

To investigate the influence of an emission-time jitter on 
the joint detection probability in a two-photon interference 
experiment, we now assume two streams of photons with a 
Gaussian emission-time distribution of identical width, Ar. 
In this case, the photon pairs are characterized by a jitter in 
the arrival-time delay of the photons, which is again given by 
a Gaussian distribution, 



1 



/n At 



cx P (-<5t 2 /At 2 ). 



(45) 



The correlation function G^(to,to + r) can be written in 
analogy to Eq. (14-11 . using only the variation of the relative 
photon delay, 

G^(t , T )= f d(6r)f(dr) trféffc, &) i(*o,*d + r)), 



(46) 

and the joint detection probability has to be calculated accord- 
ing to Eq. (I27> . In case of simultaneously impinging photons, 
this leads to 
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T 



+ At 2 



1 — cos 93 exp 



5í 2 + 5í 4 /At 2 



exp 



<5í 2 + At 2 



(47) 



In contrast to the previous case, now the width of the 
Gaussian-shaped peak in the joint detection probability of per- 
pendicular polarized photons is no more identical to the pho- 
ton duration. The variation in the emission time affects also 
the amplitude of the spatiotemporal mode functions and af- 
fects the joint detection probability even without interference. 
This can be seen in Figure[6] which shows the joint detection 
probability for photon pairs characterized by a distribution of 
the photon arrival-time delay. The width, T\, of the Gaussian- 
shaped peak is no longer identical to the photon duration, 5t. 
As one can derive from Eq. (1471 . it is now broadened by the 
width of the variation in the photon delay, Ar : 



Ti = V5t 2 + Ar 2 . 



(48) 



Furthermore, the width of the dip in case of parallel polarized 
photons is not independent of the width Ti, but is given by 



St 



T 2 = y/St 2 + <5í 4 /At 2 = — 7\. 

At 



(49) 



A time-resolved two-photon interference experiment can 
again be used to determine the variation in the emission time 
of the photon streams. However, since the shapes of the joint 
detection probabilities for a frequency and an emission-time 
jitter are identical, it is in general not possible to distinguish 
between the two. Nonetheless, one can determine the maxi- 
mum vàlues of both jitters, as well as all pairs of frequency 
and an emission-time jitters matching the data. This is dis- 
cussed in the next chapter. 

Note that in case of very short photons, the coincidence 
probability has again to be calculated using Eq.lED. The 
width of the Gaussian-shaped dip in the coincidence proba- 
bility 



cos 2 tp 



Ar 2 /ót 2 



exp 



5t 2 + At 2 



is broadened by the emission-time jitter, At, and the depth of 
the dip is also decreased by At. Again it is possible to deter- 
mine the emission-time jitter from such a two-photon interfer- 
ence experiment without time resolution. The disadvantages 
of such a procedure have already been discussed in Section 
ITVAl 



C. Autocorrelation function of the photon's shape 

We now consider only one source of single photons that we 
want to characterize using a two-photon interference experi- 
ment. We assume that the stream of photons generated by this 
source is split up in such a way that each single photon is ran- 
domly directed along two differentpaths. These two paths are 
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FIG. 6: Joint detection probability as a function of the detection time 
difference, r, for photon pairs with a variation in their relative pho- 
ton delay, St. The photons are assumed to be (a) perpendicular and 
(b) parallel polarized to each other. Curve (c) is a Gaussian mode 
function of width St. The width, Ti, of the Gaussian peak in (a) is 
broadened by the variation of the relative photon delay, Ar, which is 
here assumed to be 2 St. 



of different length and the repetition rate of the source is cho- 
sen in such a way that only successively generated photons 
impinge on the beam splitter at the same time. The details of 
such an experiment are discussed in ChapterlVl 

In general, we must distinguish between the photon ensem- 
ble of the whole stream and the subensemble of photon pairs 
superimposed on the beam splitter. The latter consists only of 
successively generated photons and the characterization of a 
single-photon source by a two-photon interference experiment 
takes into account only this subensemble. Since the jitter in 
the subensemble of successive photon pairs does not have to 
be identical to the jitter in the whole photon stream, we need a 
method to decide whether the results of a two-photon interfer- 
ence experiment can be generalized to all photons generated 
by the single-photon source. In case of very long photons, 
this can be done by comparing the joint detection probability, 
,P( 2 ) (t), for perpendicular polarized photons and the autocor- 
relation function, (t), of the average detection probability 
of the whole photon stream. 

We start our discussion with a stream of identical single 
photons, so that the average detection probability is simply 
given by the square of the amplitude of the spatiotemporal 
mode function, pW(t ) = T e 2 (t ). On the one hand, the 



12 



autocorrelation function of P^\to) reads 
AM(t) = 

fàto P«(í )P (1) (ío +r) = T 2 /dí (e(ío)e(ío +r)) 2 . 

(51) 

On the other hand, the joint detection probability for per- 
pendicular polarized photon pairs of this stream is given by 
Eq. i32i . which leads to 

P (2) M«T 2 í dí (e(í )e(í + r)) 2 . (52) 

Therefore the joint detection probability for perpendicular po- 
larized photons and the autocorrelation function have the same 
shape. 

However, if the photons show a variation in their spatiotem- 
poral modes, the two functions are no longer equal. In the 
following, we assume that the variations in the whole pho- 
ton stream are described by a normalized distribution func- 
tion /(i?), whereas the variations in the subensemble of suc- 
cessively emitted photons is given by /(#). In general, both 
functions do not have to be identical, i.e. the jitter in the 
subensemble can be smaller than the jitter in the whole photon 
stream. The autocorrelation function of the average detection 
probability P^'(to), see Eq. (I14> . is given by 

AW(t) = 

T 2 fdt JJ d^M 2 /(0i)/(0 2 ) (e(*ó, W*o + t, tf 2 )) 2 . 

(53) 

The joint detection probability, 

P( 2 )(t) oc 

T 2 Jdt JJ di?idí9 2 />i)/> 2 ) (e(í ,i?i)e(ío + T,i? 2 )) 2 

(54) 

is therefore only equal to A (2) (r), if (e(í , i?i)e(ío + t, ï? 2 )) 2 
is either independent of i?i and or if the distribution func- 
tion f(§) of the whole photon stream is identical to the distri- 
bution function /($) of successively emitted photons. 

The comparison of the joint detection probability of perpen- 
dicular polarized photon pairs and the autocorrelation func- 
tion of the average detection probability therefore answers the 
question whether the results of a two-photon interference ex- 
periment can be generalized to the whole photon stream. 

V. EXPERIMENT AND RESULTS 

In the previous three chapters we discussed the theoreti- 
cal background for characterizing single photons using two- 
photon interference. Now, we show how to use this method 
to experimentally characterize single photons that are emit- 
ted from only one source. This single-photon source has been 
realized using vacuum-stimulated Raman transitions in a sin- 
gle Rb atom located inside a high-finesse optical cavity. In 
Section lV Al we briefiy review the principle of this source and 
discuss the experimental setup, which was used to investigate 
the two-photon interference. For further details concerning 




high-finesse trigger and 

cavity recycle làser 



FIG. 7: Single-photon source and experimental setup used to inves- 
tigate the two-photon interference of successively emitted photons. 
The photons are generated in an atom-cavity system by an adiabati- 
cally driven stimulated Raman transition. A polarizing beam splitter 
directs the photons randomly into two optical fibers. The delay from 
photon to photon matches the travel-time difference in the two fibers, 
so that successively emitted photons can impinge simultaneously on 
the 50/50 beam splitter. Using a half-wave plate, the polarization of 
the photons can be chosen parallel or perpendicular to each other. 
The photons are detected using avalanche photodiodes with a detec- 
tion efficiency of 50 % and a dark-count rate of 150 Hz. 

the single-photon generation, we refer to Kuhn et al. Í22Í1 and 
references therein. The measurement of the average detection 
probability of a stream of photons emitted from this source 
is discussed in Section IV Bl As already discussed, we use 
the autocorrelation function of the average detection proba- 
bility to determine whether the results of a two-photon inter- 
ference experiment can be generalized to the whole photon 
stream. Since the duration of the photons is much longer than 
the time resolution of the detectors, the interference of succes- 
sively emitted photons is measured in a time-resolved manner. 
The results and the interpretation of these measurements are 
discussed in detail in Section lyCl and rVDl 



A. Single-photon source and experimental setup 

A sketch of the single-photon source and the experimen- 
tal setup that we use to superimpose successively gener- 
ated photons on a 50/50 beam splitter is shown in Figure 
The single-photon generation starts with 85 Rb àtoms re- 
leased from a magneto-optical trap. The àtoms enter the 
cavity mostly one-at-a-time (the probability of having more 
than one atom is negligible). Each atom is initially pre- 
pared in |e) = \5Si/2,F = 3), while the cavity is reso- 
nant with the transition between \g) = \5Si/ 2 ,F = 2) and 
\x) = \5Ps/2;F — 3). On its way through the cavity, the 
atom experiences a sequence of làser pulses that alternate be- 
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tween triggering single-photon emissions and repumping the 
atom to state |e): The 2 /is-long trigger pulses are resonant 
with the |e) <-» |af) transition and drive an adiabatic passage 
(STIRAP) to \g) by linearly increasing the Rabi frequency. 
This transition goes hand-in-hand with a photon emission. In 
the ideal case, the duration and pulse shape of each photon 
depend in a characteristic manner on the temporal shape and 
intensity of the triggering làser pulses 12411 . As we will dis- 
cuss in Section IV Cl the photon-frequency can be chosen by 
an appropriate frequency of the trigger làser 1 30] . Between 
two photon emissions, another làser pumps the atom from \g) 
to \x), from where it decays back to |e). While a single atom 
interacts with the cavity, the source generates a stream of sin- 
gle photons one-after-the-other. The efficiency of the photon 
generation is 25 %. 

As described in detail by Legero ll37ll . the source has been 
optimized with respect to jitters by compensating the Earth's 
magnètic field inside the cavity and by adding to the recycling 
scheme an additional 7r-polarized làser driving the transition 
|5<Si/2, F = 3) <-* |5F 3 / 2 , F = 2) to produce a high degree 
of spin-polarization in 55*1/2, F = 3, mp = ±3. This results 
in an increased coupling of the atom to the cavity. We have 
characterized the emitted photons by two-photon interference 
measurements before and after this optimization. 

To superimpose two successively emitted photons on the 
50/50 beam splitter, they are directed along two optical paths 
of different length. These paths are realized using two polar- 
ization maintaining optical fibers with a length of 10 m and 
1086 m, respectively. Since the photon polarization is à priori 
undefined, a polarizing beam splitter is used to direct the pho- 
tons randomly into the long or short fiber. The time between 
two trigger pulses is adjusted to match the travel-time differ- 
ence of the photons in the two fibers, which is At = 5.28 fis. 
With a probability of 25 %, two successively emitted photons 
therefore impinge on the beam splitter simultaneously. In ad- 
dition, we use a half-wave plate to adjust the mutual polariza- 
tion of the two paths. 



B. Average detection probability 

First, we investigate the average detection probability of the 
photon stream. For this measurement, the long fiber is closed 
and the detection times of about 10 3 photons are recorded with 
respect to their trigger pulses. From these photons, we cal- 
culate the probability density for a photodetection, shown in 
Figure[3] The measurement has been done (a) before and (b) 
after optimizing the single-photon source. The autocorrela- 
tion functions of both curves, shown in the inset of Figure|8] 
were calculated using Eq. i5l\ . The width of these curves is 
(a) 1.07 /ís and (b) 0.81 fis. In Section fV Dl we compare both 
results with the joint detection probability of perpendicular 
polarized photon pairs. 

Note that one obtains no information on the shape or du- 
ration of individual photons from a detection probability that 
is averaged over many photodetections. As already stated in 
Section IIV Bl such a measurement does not exclude that the 
photons are very short so that the average probability distri- 



Average detection probability 




Time [fis] 

FIG. 8: Probability density for photodetections averaged over 10 3 
photons (a) before and (b) after optimizing the single-photon source. 
The data are corrected for the detector dark-count rate. The inset 
shows the corresponding autocorrelation functions, A^ 2 \t). 

bution reflects only an emission-time jitter. Only from a time- 
resolved two-photon interference experiment, one obtains in- 
formation on the duration of the photons. This is discussed in 
the next section. 



C. Time-resolved two-photon interference 

The detection times of about 10 5 photons are registered 
by the two detectors in the output ports of the beam splitter, 
while the photons in each pair impinge simultaneously, i.e. 
with ót = 0. The number of joint photodetections, N^ 2 \ is 
then determined from the recorded detection times as a func- 
tion of the detection-time difference, t (using 48 ns to 120ns 
long time bins). To do that, the photon duration must exceed 
the time resolution of the detectors. Otherwise, joint detec- 
tion probabilities could only be examined as a function of the 
arrival-time delay, 5r, like in most other experiments. 

We have performed these two-photon interference experi- 
ments before and after optimizing the single-photon source. 
The results are shown in Figure|9]and^| respectively. Each 
experiment is first performed with photons of (a) perpendic- 
ular and then with photons of (b) parallel polarization, until 
about 10 5 photons are detected. Note that the number of joint 
photodetections is corrected for the constant background con- 
tribution stemming from detector dark counts. 

In case of perpendicular polarization, no interference takes 
place. In accordance with Section ITll FI the number of joint 
photodetections shows a Gaussian peak centered at r = 0. 
This signal is used as a reference, since any interference leads 
to a significant deviation. If we now switch to parallel polar- 
ization, identical photons are expected to leave the beam split- 
ter as a pair, so that their joint detection probability should be 
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FIG. 9: Number of joint photodetections in 120 ns-long time bins as 
a function of the detection-time difference, r, before optimizing the 
single-photon source. The results are shown for photons of (a) per- 
pendicular polarization and (b) parallel polarization. In both cases, 
the data is accumulated for a total number of 73000 photodetections. 
The sòlid lines are numerical fits of the theoretical expectations to 
the data. Ti is the width of the Gaussian peak in (a), and the width of 
the dip for parallel polarized photons (b) is given by T<x . The dotted 
curve shows the T3-wide autocorrelation function, A^ 2 '(t), of the 
average detection probability. 



FIG. 10: Number of joint photodetections in 48 ns long time bins 
after optimizing the single-photon source. The data is shown in (a) 
for perpendicular and in (b) for parallel polarized photons. In both 
cases, the data is accumulated for a total number of 139000 pho- 
todetections. The dotted curve shows the T3-wide autocorrelation 
function of the average detection probability. Compared to the re- 
sults of Figure|5| the width of the Gaussian peak in (a) is decreased 
to Ti = 0.64 /iS and it is clearly smaller than T 3 = 0.81 fj,s. With 
parallel polarizations, a dip of increased width, T2 = 0.44 us, is 
observed. All results are discussed in Section lV Dl 



zero for all vàlues of r. In the experiment, however, the signal 
does not vanish completely. Instead, we observe a pronounced 
minimum around r = 0, which complies well with the behav- 
ior one expects for varying spatio-temporal modes, as shown 
in Figures [5] and [6] In analogy to Eq. J42i and A47L respec- 
tively, the number of joint photodetections is then given by 



= 4 2) 



1 - cos tp exp ( 



(2) 

Nq ' is the peak value at r 
ular polarized photons, tp - 



(55) 

= that we measure for perpendic- 
7r/2. The time T\ is the width of 



this Gaussian peak, whereas T2 gives the width of the dip for 
photons of parallel polarization, ip = 0. In Section lV Dl both 
numbers are used to deduce the frequency and the emission- 
time jitter. The two widths, T\ and T%, are obtained from a 

fit of Eq. ( I55> to the measured data. This is done in two steps. 

(2) 

First, we obtain T\ and Nq from a fit to the data taken with 
perpendicular polarized photons. We then keep these two vàl- 
ues and obtain the dip width T2 from a subsequent fit to the 
data with parallel polarized photons. In this second step, we 
use a polarization term of cos 2 tp = 0.92 to take into account 
that we have a small geomètric mode mismatch. This is well 
justified since mode mismatch and non-parallel polarizations 
affect the signal in the same manner. The value of tp has been 
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FIG. 1 1 : Two-photon interference of parallel polarized photons with 
a frequency difference, A, of (a) 2ir x 2.8 MHz and (b) 2tt x 3.8 MHz. 
The number of joint photodetections is accumulated over (a) 210000 
and (b) 319000 detection events. It oscillates as a function of the 
detection-time difference. The sòlid curves represent numerical fits 
to the data with a frequency-difference of (a) A = 2n x 2.86 MHz 
and (b) A = 27T x 3.66 MHz. The dotted curve shows the reference 
signal measured with perpendicular polarized photons. 



obtained from an independent second-order interference mea- 
surement. 

As one can see by comparing Figures |9] and [KÏ] the com- 
pensation of the Earth's magnètic field and the improved re- 
cycling scheme lead to a decreased width Ti of the Gaussian 
peak in (a) and a broader dip T2 in (b). As we discuss in the 
following, these results show that this optimization has suc- 
cessfully reduced the jitter in the mode function of the pho- 
tons. 

Moreover, as shown in Figure II ll a) and (b), we resolve a 
pronounced oscillation in the number of joint photodetections 
when a frequency difference, A, is deliberately introduced be- 
tween the interfering photons 1 30] . This is achieved by driv- 



ing the atom-cavity system with a sequence of trigger pulses 
that alternate between two freqüències. The frequency differ- 
ence between consecutive pulses is either (a) 2tt x 2.8 MHz 
or (b) 2-7T x 3.8 MHz. In accordance with Section IIÏTFl the 
oscillation in the joint detection probability always starts with 
a minimum at r = 0, and the màxima exceed the reference 
signal that we measure with perpendicular polarized photons. 

The latter experiment has been performed with the opti- 
mized single-photon source. The number of detected photons 
(corrected for the number of dark counts) equals the photon 
number in Figure ^0| Therefore Nq 2 \Ti and T2 are well 
known from this previous measurement. As shown by Leg- 
ero 13711 . the only remaining parameter one can obtain from a 
fit of the joint detection probability, 



A^ 2 )(t) = 

(2) _ 



N, 



exp 



1 — cos 2 ip cos (t A) cxp 

(56) 

to the data is the frequency difference A. Its fit value agrees 
very well with the frequency differences that we imposed on 
consecutive pulses. We therefore conclude that the adiabatic 
Raman transition we use to generate the photons allows us 
to adjust the single-photon frequency. Moreover, the oscil- 
lations in the joint detection probability impressively demon- 
strate that time-resolved two-photon interference experiments 
are able to reveal small phase variations between the photons, 
like, e.g., the frequency difference we have deliberately im- 
posed here. 



D. Interpretation of the results 

We start our analysis by comparing the autocorrelation 
function of the average detection probability with the result 
of the two-photon coincidence measurement with perpendic- 
ular polarized photons, shown in Figures[9ja) and llOf a). The 
shape of the autocorrelation function, (r), is commensu- 
rable with the data obtained in the two-photon correlation ex- 
periment before the optimization of the source, but it differs 
significantly afterwards. As discussed in Section lïV Cl the dif- 
ferent widths of both curves, T\ = 0.64 /ís and T3 = 0.81 fis, 
indicate that the photon stream is subject to variations in the 
spatiotemporal mode functions that are much less pronounced 
in the subensemble of consecutive photons. These variations 
cannot be attributed to a jitter in the photon frequency, since 
the autocorrelation function and the joint detection probabil- 
ity, given by Eq. d53l and i54\ , depend only on the frequency- 
independent amplitude of the mode functions. Therefore the 
emission time and/or the duration of the photons must be sub- 
ject to a jitter. As a consequence, the average detection prob- 
ability shown in Figure [8] cannot represent the shape of the 
underlying single-photon wavepackets. In particular the width 
of the measured photon detection probability is broadened due 
to the emission-time jitter. 

Moreover, the discrepancy between (r) (Figure[8} and 
the Gaussian peak (Figures [9] (a) and II 01 (al) shows that the 
variations in the whole photon stream are larger than the vari- 
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FIG. 12: Frequency jitter, 5u>, and emission-time jitter, Ar, before 
(a) and after (b) optimizing the single-photon source. The two curves 
represent all pairs of jitters that match the two widths, Ti and T2, 
found in the two-photon interference experiments. 



Ar = 0.82 /zs to Ar = 0.53 /zs. 

However, in general, both the frequency and the emission 
time are subject to a jitter. If these fiuctuations are uncorre- 
lated, a whole range of (Süj, Ar)-pairs can explain the peak 
and dip widths, T\ and T%. For our two sets of data (before (a) 
and after (b) the optimization of the source), this is illustrated 
in FigurefT2l All pairs of frequency and emission-time jitters 
that are in agreement with the measured vàlues of T\ and T2 
lie on one of the two sòlid lines. From this figure, it is evident 
that the vàlues for Soj and Ar deduced above represent the up- 
per limits for the respective fiuctuations. Moreover, it is also 
nicely visible that our optimization of the source significantly 
improved the frequency stability and emission-time accuracy 
of our single-photon source. 

We emphasize again that this information about the pho- 
tons can be obtained from time-resolved two-photon interfer- 
ence experiments, but not from a measurement of the average 
detection probability. 



VI. CONCLUSION 



ations in the subensemble of consecutive photons. Therefore 
the following analysis of the times T\ and T2 of the two- 
photon interference cannot be generalized to the whole photon 
stream. 

To figure out the photon characteristics that can explain the 
measured quantum-beat signal, we restrict the analysis of Ti 
and T2 to a frequency and an emission-time jitter. If we as- 
sume that successively emitted photons show only a variation 
of their freqüències, then Tj is identical to the photon du- 
ration, whereas T2 = 2 /Su> is solely due to the frequency 
variation. In this case, T2 is identical to the coherence time, 
which one could also measure using second-order interference 
1271 12811 . As shown in Figuresl9land llOI optimizing the single- 
photon source reduces the bandwidth of the frequency varia- 
tion from Soj/2-k = 1.03 MHz to 6w/2n = 720kHz. The re- 
maining inhomogeneous broadening of the photon frequency 
can be attributed to several technical reasons. First, static and 
fluctuating magnètic fields affect the energies of the Zeeman 
sublevels and spread the photon freqüències over a range of 
160 kHz. Second the trigger làser has a linewidth of 50 kHz, 
which is also mapped to the photons. And third, diabatically 
generated photons lead to an additional broadening. 

Another explanation for the measured quantum-beat signal 
assumes photons of fixed frequency and shape, but with an 
emission-time jitter. In this case, we have T\ = \jSt 2 + At 2 
and T2 = T\6t/ At. From these two equations, one can cal- 
culate St, which is the lower limit of the photon duration, and 
Ar, which is the maximum emission-time jitter. In our ex- 
periment, the optimization of the single-photon source led to 
an increase of St from 0.29 /zs to 0.36 /zs, and at the same 
time to a reduction of the maximum emission-time jitter from 



We have shown that time-resolved two-photon interference 
experiments are an excellent tool to characterize single pho- 
tons. In these experiments, two photons are superimposed on 
a beam splitter and the joint detection probability in the two 
output ports of the beam splitter is measured as a function of 
the detection-time difference of the photons. This is only pos- 
sible if the photons are long compared to the detector time 
resolution. For identical photons, the joint detection proba- 
bility is expected to be zero. Variations of the spatiotempo- 
ral modes of the photons lead to joint photodetections except 
for zero detection-time difference. Therefore the joint detec- 
tion probability shows a pronounced dip. From the width of 
this dip, one can estimate the maximum emission-time jitter 
and the minimum coherence time of the photons. In addition, 
a lower limit of the single-photon duration can be obtained. 
This is not possible by just measuring the average detection 
probability with respect to the trigger producing the photons. 
Moreover, we have shown that a frequency difference between 
photons leads to a distinct oscillation in the joint detection 
probability. This does not only demonstrate that we are able 
to adjust the freqüències of the photons emitted from a single- 
photon source, but also that one is sensitive to very small fre- 
quency differences in time-resolved two-photon interference 
measurements. 
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